function [post,like,lik1,lik2] = fnmhsv(para,data,rho,sigu,sigrho);
global msel datasel     

     lpri =  logprior(para); 
     like =  logliksv(para,data,rho,sigu); 
     lik1 = 0;
     lik2 = 0;     
     siguw   = para(5); 
     sigrhow = para(6); 
     nv      = length(sigu); 
     
     lik1 = lik1-0.5*(log(2*pi)+log(siguw^2))-(log(sigu(1))-para(3))^2/(2*siguw^2);  
     lik2 = lik2-0.5*(log(2*pi)+log(sigrhow^2))-(log(sigrho(1))-para(4))^2/(2*sigrhow^2);  

        for t=2:nv 
     lik1 = lik1 -0.5*(log(2*pi)+log(siguw^2))-(log(sigu(t))-log(sigu(t-1)))^2/(2*siguw^2); 
     lik2 = lik2 -0.5*(log(2*pi)+log(sigrhow^2))-(log(sigrho(t))-log(sigrho(t-1)))^2/(2*sigrhow^2);   
        end 
     
   
     post  = lpri+like+lik1+lik2; 
     
